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Abstract 

It has been controversial whether rjb can be discovered in Tevatron Run 2 through the decay 
Vb J/ip J/ip followed by J/tp — > u + u~ . We clear this controversy by an explicit calculation which 
predicts Br [776 — ► J/ifi J/ifi] to be of order 10 -8 . It is concluded that observing rjb through this decay 
mode in Tevatron Run 2 may be rather unrealistic. The rjb may be observed in the forthcoming 
LHC experiments through the 4-lepton channel, if the background events can be significantly 
reduced by imposing some kinematical cuts. By some rough but plausible considerations, we find 
that the analogous decay processes rjb — > VV, D*D also have very suppressed branching ratios, 
nevertheless it may be worth looking for rjb at LHC and Super B factory through the decay modes 
rib -» KsK^ir?, D*D. 
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I. INTRODUCTION 



Since the discovery of T(1S) nearly three decades ago, extensive search for its pseudo- 
scalar partner, rjb, has been conducted in various experiments. Unfortunately, to date there 
is still no conclusive evidence that this elusive particle has been found [lj]. 

The existence of rjb is a solid prediction of QCD. On the theoretical side, many work have 
attempted to unravel its various properties. In particular, its mass is believed to be among 
the simplest and most tractable observables. Numerous estimates for T — 775 mass splitting 
span the range 20-140 MeV , [h], [ll|. Among different theoretical 

approaches, perturbative QCD is believed by many people to yield reliable predictions, 
because of decently heavy b mass and rjb being the lowest-lying bb state. By far the most 
sophisticated prediction along this direction, facilitated by the NRQCD renormalization 
group technique, gives M Vb = 9.421 ±0.013 GeV [h]]. An eventual unambiguous observation 
of rjb and precise measurement of its mass will decisively test the weakly-coupled picture of 
the bb ground state. 

Much efforts are spent to search rjb from 77 collisions in the full data samples of LEP 
2, where approximately two hundreds of rjb are expected to be produced. ALEPH has one 
candidate event rjb — > -Ks(— > 7t + 7t~)K~7t + 7t~tt + (possibly missing a 7r°) with the recon- 



structed mass of 9.30 ± 0.03 GeV, but consistent to be a background event 



12| . ALEPH, 



L3, DELPHI have also set upper limits on the branching fractions for the rjb decays into 4, 



6, 8 charged particles 
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141 ]. Based on the 2.4 fb 1 data taken at the T(2S) and T(3S) 



resonances, CLEO has searched distinctive single photons from hindered Ml transitions 
T(2S), T(3S f ) — > r7;,7 and T(3S) — > hbTt , hbn + n^ followed by El transition hb — > TjbJ, but 
no signals have been seen }15| . 

Hadron collider experiments provide an alternative environment to search for rjb- Unlike 
the e + e~ machines which is limited by the low yield of rjb, hadron colliders generally possess 
a much larger production cross sections for 776, which in turn allows for searching it through 
some relatively rare decay modes yet with clean signature. However, one should bear in 
mind that a noteworthy disadvantage also accompanies with hadron collision experiment, 
i.e., that the corresponding background events may also be copious, so the effectiveness of 
these decay modes might be seriously discounted (Such an example may be the decay mode 
T)b —> 77, with branching fraction ~ 10~ 4 , but the combinatorial background 7 events can 
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be enormous). 

Encouraged by the large observed width of r) c — ► VV(V stands for light vector mesons), 
Braaten, Fleming and Leibovich (hereafter BFL) have suggested that the analogous decay 
process rjb — ► J/ipJ/ip, followed by both J ftp decays to muon pairs, may be used as a very 



clean trigger to search for rjb at Tevatron Run 2 16j. Assuming Br [rjb — > J/ipJ/ip] ~ 

they rescale the measured branching ratio of rj c — > 00 by a factor of (m c /m;,) 4 to estimate 1 

Br[% -> J/ipJ/ip] = 7 x 10" 4±1 . (1) 

Combining with the knowledge about the production rate of rjb at Tevatron, BFL conclude 

that the prospect of observing rjb through the 4/x decay mode in Run 2 is promising. 

Following this suggestion, CDF has searched for the rjb — ► J/ipJ/ip — > 4/x events in the 
I I 

full Run 1 data sample [19j]. A small cluster of seven events are seen in the search window, 
where 1.8 events are expected from background, with the statistical significance of 2.2 a. A 
simple fit infers the cluster's mass to be 9.445 ±0.006 (stat) GeV. If this cluster is truly due 
to r7b, then the product of its production cross section and decay branching ratio is close to 
the upper end of BFL expectation. 

In a recent work, Maltoni and Polosa (henceforth MP) nevertheless argue that the BFL 
estimate, (pQ), may be overly optimistic 20j. MP suspect that the analogy between rj c — > 00 
and rjb — > J/ipJ/ijj is only superficial. It is known that perturbative QCD (pQCD) framework 
has difficulty to account for the large observed widths of the r] c VV decay processes |2ll . 



red to reconcile 
23, hi. On the 



221 ] . The consensus is that some nonperturbative mechanisms should be invo 
the discrepancy between the pQCD prediction and the actual measurement 
other hand, due to heavy b and c masses, it is rather reasonable to expect that rjb — > J/ipJ/ip 
can be safely tackled within the pQCD scheme. Therefore, the rescaling procedure used by 
BFL, whose validity should reside only in the domain of pQCD, may well be illegitimate 
when taking the measured ratio of r) c — > 00 as input, which is essentially dictated by 
nonperturbative dynamics. 

The persuasive evidence in favor of MP's argument comes from their explicit calculation 
for the inclusive decay rate of rjb to 4-charm states, 

Bi[rj b -> cccc] = l.StH x 10~ 5 . (2) 



1 Note that Br GX p[?y c -> <j><j>] has shifted from (7.1 ± 2.8) x 1CT 3 given in 2000 PDG edition [lJJ, which was 
quoted by BFL, to the latest PDG value (2.7 ± 0.9) x 10" 3 [18 1. 
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This ratio is even smaller than the lower limit of Br [77;, — > J/ipJ/ip] estimated by BFL! 
This is clearly at odds with the usual thought that the exclusive decay rate should be much 
smaller than the inclusive one. 

One of the major concerns of this paper is to dispose of this controversy by performing an 
explicit calculation. Heavy b and c quark masses set hard scales so that we can confidently 
utilize pQCD to tackle this decay process, expecting those nonperturbative contributions 
plaguing the decay r] c — > VV play an insignificant role here. Since each involved particle is 



a heavy quarkonium, we work with the color-singlet model, in line with t 



26 



27 



28 



re calculation done 



29 



30 



31| . It is found 



for double charmonium production at e + e~ colliders 
that at the lowest order in a s , retaining the transverse momentum of c inside J/ip is vital 
to obtain a non- vanishing result, and consequently, the correct asymptotic behavior of the 
hadronic decay branching ratio is a^vl°(m c /mb) 8 (v c stands for the typical velocity of c in 
J/ip). Numerically, we predict 



Br[?7 6 -»■ J/ipJ/ip\ = (0.5 - 6.6) x 10~ 8 , (3) 

which is much smaller than the BFL estimate, ([1]). Simple analysis indicates that the cluster 
reported by CDF [l9] is extremely unlikely to be affiliated with rjb- We further argue that, 
the potential of discovering rjf, through this decay mode is gloomy even in Run 2. 

To better guide experimental search for rjb, an important issue is to know that, besides 
Vb — > J/^J/^y which hadronic decay modes can serve as efficient triggers to detect 
Heavy 77^ mass opens a huge phase space available for innumerable decay modes, but at a 
price that the branching ratio of each individual mode has been greatly diluted with respect 
to that from rj c decay. 

At any rate, it is valuable to know which hadronic decay modes of 775 possess the largest 
branching ratios. We examine numerous hadronic decay channels, mostly two-body ones, 
e.g., r] b decays to two light mesons, and two charmed mesons. By some crude but plausible 
estimate, we find 77^, — > 00, D*D* also have very tiny branching ratios, of the same order as 
that for rjf) — > J/ip J/ip. In sharp contrast, the branching ratios for 775 — > D*D might be as 
large as 10 -5 , therefore they can be used as searching modes. Furthermore, stimulated by 
the experimental fact that rj c decays to three pseudoscalars have largest branching ratios, 
we urge experimentalists to look into such 3-body channels as 77^ — > KKir, 77(7/')^^, DDtt. 
The corresponding ratios are estimated to be of order 10~ 4 . 
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FIG. 1: Lowest-order QCD diagrams which contribute to 775 — > J/tp J /if). 



The remainder of the paper will proceed as follows. In Sec. [TTl we discuss the asymptotic 
behavior of Br [rji, — ► J /if) J /if)] in the limit 7715 ^> m c ^> Aqcd ; and elucidate the peculiarity of 
this decay process. We then present the actual calculation of the ratio in Sec. IHIl employing 
the color-singlet model which incorporates velocity expansion. In Sec. IIVI we discuss the 
observation potential of rjb through this mode in Tevatron and LHC In Sec. |V] by some 
simple scaling analysis, we estimate the branching ratios for rj^ decays to various final states, 
e.g., to two light mesons, to three light pseudo-scalar mesons, and to two charmed mesons. 
For the decay process r/?, — > VV, we also compare our estimates with those obtained from 
some nonperturbative mechanism. In Sec. IVII we summarize and give an outlook. 



II. ASYMPTOTIC BEHAVIOR AND UNNATURAL DECAY PROCESS 



Before launching into the actual calculation, it is instructive to first envisage the general 
property of this exclusive decay process. A powerful tool in pQCD to count kinematical 
suppression factors for hard exclusive processes is the so-called hadron helicity selection 
rule, originally developed for the reactions involving light hadrons [32]. This rule has re- 
cently been applied to analyze the double charmonium production in e + e~ annihilation in 



the limit y/s ^> m c 25j. One can work out the asymptotic behavior for Br[^ — > J /if) J /if)] 



in an analogous way. The decay 775 — > J /if) J /if) can be initiated by either strong or elec- 
tromagnetic interaction, with the corresponding lowest-order diagrams shown in Fig. [1] and 
Fig. [2, respectively. Generally speaking, the QCD contribution is dominating over the QED 
contribution. However, for the sake of completeness, the latter will be also included in our 
analysis. 
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For the lowest order strong decay process depicted in Fig. [TJ a simple consideration may 
suggest that in the limit m& — > oo with m c fixed, 

/m 2 \ 2+|A+X| 

Br str [ % - J/^(A) + J/V(A)] ~ a2 «2 -| , (4) 

where A and A represent the helicities of two J ftp viewed in the rji, rest frame, and v c 
denotes the characteristic velocity of c in J/ip. Obviously the scaling behavior depends 
on the helicity configurations of both J/i/j. The factor v\ is expected because there is a 
factor of wave function at the origin, ^ju(O), for each cc pair to emerge with small relative 
momentum to form S-wave bound state, and if)j/^(Q) ~ (m c t> c ) 3//2 . 

The expectation (j3J) is compatible with the helicity selection rule that the decay config- 
uration which conserves the hadron helicity, i.e. A + A = 0, exhibits the slowest asymptotic 
decrease, Br str ~ l/ m &- This occurs because there are two hard gluons carrying large mo- 
mentum transfers. The only helicity state bearing this least suppressed ratio which is also 
compatible with the angular momentum conservation A = A is thus (A, A) = (0, 0). The he- 
licity conservation can be violated either by the nonzero charm mass m c or by the transverse 
momentum of c inside J/ip, q±. For every unit of violation of the selection rule, there is 
a further suppression factor of m 2 c jm\ or q^/m 2 .. For the other physically allowed configu- 
rations (A, A) = (±l, ±1), the helicity conservation is violated by two units, so one expects 
Br str ~ 1/mf. By writing (j3J) the way as it is, we have temporarily assumed that the cause 
of violation is entirely due to the quark mass m c . 

The essential assumption of BFL is the "leading twist" scaling behavior Br str ~ 
which is tacitly associated with the T/& decay to two longitudinally polarized J/ip 2 . However, 
one may recall that rjf, — > J/ifjJ/ip, like r\ c — > VV, belongs to a class of so-called unnatural 
decay processes |33(, for which the helicity state (0, 0) is strictly forbidden due to the conflict 
between parity and angular momentum conservation 3 . In operational basis, it arises because 
the decay amplitude for such a process involves the Levi-Civita tensor and there are no 
enough number of independent Lorentz vectors to contract with it for vector mesons are 



2 The decay r\h — > J/ipr) c would exhibit the "leading twist" scaling behavior, Br ~ o? s v^(m c /mb)^, were it 
not inhibited by C invariance. 

3 An unnatural process, according to Ref. 33|, is defined as a heavy quarkonium two-meson decay process 



that does not conserve a multiplicative quantum number called naturalness, which is defined by a 
(— 1) J P (J, P stand for the spin and parity of a meson). 
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longitudinally polarized. Further examples of the unnatural processes for a bottomonium 
decays to two S- wave charmonia include T {xhst) J/^Vc and Xbi J/ipJ/ip, in contrast 



with the natural decay processes such as Xbo,2 — > J/ipJ/ip |34| . r) c r] c , and hb — > J/ipr] c . 

We may wish to examine this assertion closely for our case. Parity and Lorentz invariance 
constrain the decay amplitude to have the following tensor structure: 

M(X,X) oc e^P M P,4(A)^(A), (5) 

where P, P, e and e are momenta and polarization vectors for both J/ip. If both J/ip are 
longitudinally polarized, e and e then can be expressed as linear combinations of P and P, 
M(0, 0) thus vanishes 4 . It is worth emphasizing that this result is based solely on the basic 
principle of parity and angular momentum conservation, so it will not depend on dynamical 
details. For example, it will be true irrespective of whether this process is initiated by QCD 
or QED, whether higher order perturbative corrections are included or not, and whether 
nonperturbative QCD effects are incorporated or not. 

In passing, one may notice an equivalent but more intuitive explanation for M(0, 0) = 



for the unnatural decay processes [22J. It can be attributed to a peculiar property of Clebsch- 
Gordon coefficients that (10|10; 10) = 0. We take rjb — > J/tpJ/ip as an example to illustrate 
this. Due to parity and angular momentum conservation, two J/ip must have the relative 
orbital angular momentum L = 1 and the total spin 5 = 1. However, it is impossible for 
the two longitudinally polarized J/ip to couple to a (S, S z ) = (1, 0) state because of this 
vanishing Clebsch-Gordon coefficient. 

Since both J/ip in the final state must be transversely polarized, the helicity conservation 
is violated by two units, Eq. (J4]) then indicates Br str [?7b — > J/4>± J/ip±] ~ a 2 s v® (m c /mb) 8 . 
Nevertheless, as the explicit calculation will reveal, this behavior is not quite correct in 
counting powers of v c . It turns out that the true asymptotic behavior is even more sup- 
pressed, 

Br stI [ Vb - m±l) + m±l)] - * 2 A ~ alvf (^) 8 , (6) 

where q± ~ m c v c is assumed in the last term. This implies that, at the lowest order in a s , the 
violation of the rule should be ascribable to the nonzero transverse momentum of c in J/ip, 



4 Exactly with the same argument, the decay e + e — > 7* — > J/tp(X = 0) + ?/ c turns out to be strictly 
forbidden. 
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FIG. 2: Lowest-order QED diagrams that contributes to —* J/ip J/tp- Only the fragmentation 
type diagrams are retained, whereas the other two, which can be obtained by replacing the gluons 
in Fig. [TJby photons, have been suppressed. 



instead of the nonzero m c . This is compatible with the earlier finding that the amplitude 
for f] c — > V± V± vanishes in the collinear quark configuration, even though the quark masses 



are kept nonzero 



211 ] . After q±_ is included for the light vector mesons, the decay rate then 



becomes nonzero, though still too tiny to account for the measured rates 

We can also infer the asymptotic behavior of the electromagnetic contribution in Fig. [21 
The photon fragmentation produces transversely polarized J/ip. Since the hard scale is set 
by the virtuality of fragmenting photon ~ m 2 cl no suppression factor oc l/mf can arise in the 
branching ratio. A simple counting rule suggests that the QED fragmentation contribution 
to the ratio exhibits the following behavior 

Br em [ % - JM±1) + JM±1)] ~ (-^ vt . (7) 



Although the electromagnetic contribution is free of suppression by inverse powers of m^, in 
general it still cannot counteract the adversity caused by a < a s . 

There is also the interference term between QCD and QED contributions, with a scaling 
behavior intermediate between (ED and (TTl). 



III. COLOR-SINGLET MODEL CALCULATION 

In this section, we present a calculation for rjb — > J/i/jJ/i/j in perturbative QCD scheme. 
As stressed before, this scheme is expected to generate reliable answer for this process, since 
the annihilation of bb and creations of cc pairs take place in rather short distances. This 
decay process should be contrasted with the analogous process r) c — > VV. To account for the 
unnaturally large measured width of the latter, it is proposed that some nonperturbative 
mechanisms, e.g., the mixing among pseudo-scalar states rj — rj' — rj C) either due to QCD axial 
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anomaly [23j, or due to perturbative box diagram [24(, together with quark pair creation from 
vacuum, should play a prominent role. The impact of final state interaction on the doubly- 
OZI-suppressed process r] c — > uj<p has also been addressed j^]. Fortunately, neither of these 
nagging complications will bother us, because the //& is too heavy to bear a significant mixing 
with t] c , and c is too heavy to easily pop out of the vacuum. The weak inter-quarkonium 



Van der Vaals interaction 36] implies that the final state effects are also unimportant. 

We first consider the QCD contribution in Fig. [I] for r)b(K) — > J/ip(P) + J/ip(P), where 
K, P and P signify momenta of each quarkonium. In this work, we will only consider 
the color-singlet Fock state of the quarkonium, completely ignoring the possible color-octet 
effect, which is difficult to analyze for exclusive processes in a clear-cut way 5 . 

In the color-singlet model calculation, it is customary to begin with the parton level 
matrix element b(k)b(k) — > c(p)c(p) + c(p)c(p), then project this matrix element onto the 
corresponding color-singlet quarkonium Fock states. It is worth stressing that because of 
impossibility to find a frame such that two back-to-back fast moving J/ip become simulta- 
neously at rest, it is imperative to admit a manifestly Lorentz-covariant projector. 

For the cc pair with total momentum P, we assign the momentum carried by each con- 
stituent as 

P 

V = 77 + Q, 
P 

P = 2" -?> ( 8 ) 

where q is the relative momentum satisfying P ■ q = 0. In the rest frame of cc pair, P^ is 
thus purely time-like and q^ purely space-like. Since the cc pair forms J/ip, it is necessarily 
in a spin-triplet color-singlet state, and one can replace the product of the Dirac and color 



spinors for c and c in the final state with the following Lorentz-covariant projector 
v{p)u{p) — |[d ff ]—^—l_ -tf-m) f(X)(P + 2E)(^+m) 



25 



38|: 



x 



1 , , o.\ ( 1 



w (-g 2 ) ® . (9) 



m ' \ JN, 



Note the color-singlet model can be viewed as a truncated version of the Nonrelativistic QCD (NRQCD) 



factorization approach 37|. which still admits a factorization of the calculable short-distance part and 
universal long-distance factors in the color-singlet channel. For the exclusive processes involving heavy 
quakonium, the color-single model and NRQCD are often used interchangeably in literature. 



9 



where p 2 =p 2 = m 2 , E = yP 2 /2 = Jm 2 — q 2 is a Lorentz scalar, and e* is the polarization 
vector of J/ip satisfying e(X) • e*(A') = — 5 XX ' and P ■ e* = 0. N c = 3 and l c stands for the 
unit color matrix. The momentum space wave function (pj/^pi—q 2 ) is explicitly included and 
the Lorentz-invariant measure is defined as 

The introduction of the covariant projector (jHj) unites the non-relativistic internal motion 
and highly relativistic external motion of J/ip in a coherent fashion, which is indispensable if 
one wishes to systematically implement relativistic corrections for highly energetic processes 
like this one. At the final stage of the calculation, one can always choose to perform the 
momentum integrals in the rest frame of each J/if> successively, thanks to the Lorentz- 
invariant measure (110p . We take the following integral, which arises in the zeroth order of 
relativistic expansion, as an example: 



J[dq](pj/i>(-q 2 ) = Jj^^J/i*^ 



= V\w(o), (11) 

rest frame 



where i[>j/$(0) is the spatial Schrodinger wave function at the origin for J/ip. 

The above formulas obviously also apply to the second J/ip, once the replacements P — ► 
P, q — > q, E — > E, and e(A) — > e(A) are made. 

For the 775 in the initial state, it turns out that the 0(v 2 ) correction is less relevant, so 
we will neglect the relative momentum and simply take k = k = K/2. Consequently, the 
following simplified projector will be used: 

u(k)v(k) — ► -L(?( + 2m b )z l5 x ^^^(0)) ® (4=) ■ (12) 



c / 

For processes involving heavy quarkonium, one usually organizes the amplitude in powers 
of the relative momenta, to accommodate the NRQCD ansatz. Some simple algebra shows 
that the rjb part of the matrix element contains the factor 

Tr[ 75 W(4- fhl = Ue^K^q-q)*. (13) 

In the leading order in v c expansion, q = q = 0, the amplitude hence vanishes. As a matter 
of fact, Anselmino et al have drew the same conclusion for r] c — > VV while using light-cone 



scheme in the collinear quark configuration but retaining nonzero light quark masses 



2i| 
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This color-singlet-model result can be viewed as a specific example of theirs by invoking 
narrow-peak approximation to the light-cone wave function for vector mesons. 

In order to obtain a non-vanishing amplitude, one must go to the next-to-leading order 
in v c . In particular, one should expand the remaining part of amplitude to linear order in q 
and q, to pair up with the ones in Eq. flTSJ) 6 . 

First we expand the product of three propagators in Fig. [TJ 

1 1 1 



(q - q) 2 - m 2 {K/2 - q + q) 2 (K/2 + q - q) 2 



1 



rn, 



1 + 



q 2 q ■ q 



(f 



1 ' 2 ' 2 

mi m^ m^ 



(14) 



Because the sub-leading terms are at least quadratic in q or q, they can be dropped. 

The missing q or q thus must arise from the J/ip part of the matrix element, or more 
precisely, from the projectors for J/ip, Eq. (Q. After some straightforward Dirac trace 
algebra, one can express the full QCD amplitude in the Lorentz-invariant form: 



M 



str 



m 



13/2 



(15) 



where the color factor C str 
integral A reads 

A = 

x 



iV- 3 / 2 Tr (T a T b ) Ti (T a T b ) = A^ 1 N~ 3 ^ 2 (N 2 - 1), and the double 



,~2\ 



[dq] [dq\(j>j/^{-q )(j)j/^(-q 



e^PvP" q a {a*? q-E*-e*Pq- £*) + (q q) 



(16) 



where only the terms with the intended accuracy of 0(v 2 ) are kept. We have substituted 
E, E everywhere by m c , which is legitimate in the current level of accuracy 7 . 

One immediately observes that, to survive when contracted with the antisymmetric ten- 
sor, q and q, as well as e* and e* in (fl6|) . must be transverse. We thus confirm the earlier 
assertion that transverse momentum is the agent to violate the helicity conservation. 

Because the transverse components of a 4- vector are invariant under the boost along the 
moving direction of J/if), we can perform the integrals in the rest frame of J/ip without 



6 Retaining the relative momentum in 775 but not q and q still leads to vanishing amplitude. 

7 For the same reason, we have not bothered to include explicitly the relativistic effects due to 2-body phase 
space and normalization of cc states as considered in Ref. [2d. |38| . which are powers of E/m c or E/m c . 
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concerning about the boost effect. Using the spherical symmetry of S'-wave wave function, 
we have 



f f d?q 

J[dq}q t ± q{<i ) J/4- ( l 2 ) = y J J^y ^ 

ft 



rest frame 



= y^{« 2 )j/^ W (0). (17) 

Here (v 2 ) jm is a quantity governing the size of relativistic corrections. Loosely speaking, it is 
related to the second derivative of wave function at the origin for J/ip, and characterizes the 
average v 2 of c inside J/ip- Inspecting Eq. ([17]) . however, one immediately realizes at large 
q, the Coulomb wave function should dominate and this quantity turns out to be linearly 
ultraviolet divergent, hence its meaning becomes obscure. In fact, this factor admits a more 



rigorous definition as a ratio of NRQCD matrix elements 25|, [371] : 



, 2 , e-(J/^e)\^a(-B 2 ) X \0) 

{ )m mle-(me)W<T X \0) ' ^ 
where ip and x represent Pauli spinor fields in NRQCD, D is spatial covariant derivative. The 
matrix elements appearing in the above ratio should be understood to be the renormalized 
ones. Lattice QCD extraction of this quantity has been available long ago, but the precision 



391 ] . We will specify our choice for the numerical value of (v 2 )j/^ in Section 



is quite poor 

ED 

Substituting (|T7|) into f|T6|) . we obtain 



A = -iml{f) J/ 4ift / i(Q)e^P*P»e* a i* , (19) 
which has the desired tensor structure as in (jHJ). The strong decay amplitude then reads 

M stt = 512 ^*]$ m ° ^ b (o) V^(Q) {v 2 )m p* p* < n ■ (20) 

27 m b 

Next we turn to the electromagnetic contribution to rjb — > J/ipJ/ip. Two QED diagrams 
which have the same topology as Fig. [fl but with gluons replaced by photons, lead to 
the amplitude of the same form as (12"U|) except a 2 is replaced by el e 2 c a 2 . Obviously, their 
contributions are much more suppressed than those from the fragmentation diagrams in 
Fig. |2l hence will not be considered. The QED fragmentation contribution to the amplitude 
can be easily worked out, 

24 ^6 7T 2 e 2 e 2 a 2 , _ , 9 _ 2m 2 

5/2 r. 



M em = " V 7 /2 e 7 a ^(0)^(0) 1-^ e^P.P^. (21) 

m 6 \ Tfli, / 
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Adding (|20|) and (I2TT) together, squaring, summing over transverse polarizations of both 
J/ip and integrating over half of the phase space, we then obtain the partial width — > 
J/ipJ/ip}. Nevertheless, it is more convenient to have a direct expression for the branching 
ratio, where ip Vb {0) drops out, 



Br[r/ ft -> J/iP J/iP] 



2 13 7i 2 a 2 m 



x 



3 4 

,2 



1 - 



4m 



2\ 3/2 



m 



m 



'9 e& e c a m| x 2 



•a. 



ml? J \ 



2m;? 
m 2 



In deriving this, we have approximated the total width of rji, by its gluonic width: 



57T<X 



3 m 2 



(22) 



(23) 



where the LO expression in a s and v b is used for simplicity. 

Equation (|22p constitutes the main formula of this paper. One can readily confirm the 
asymptotic behavior of QCD and QED contributions first given in Sec. [TTJ Eq. fl6]) and (J7|). 



IV. OBSERVATION POTENTIAL OF rj b -> J/ip J/ip AT TEVATRON AND LHC 



We now explore the phenomenological implication of (1221) . The input parameters are mj, 
m c , a, a s , ipj/^(0) and {v 2 ) j/^, all of which can be inferred from other independent sources. 
The wave function at the origin for J/ip can be extracted from its dielectron width: 

, , 47re 2 a 2 



r[J/v 



e e 



mt 



(24) 



(The LO formula in a s and v 2 is used for simplicity.) Using the measured dielectron width 
5.55 keV 3, we obtain ipj/^O) = 0.205 GeV 3/2 for m c = 1.5 GeV. 

Among various input parameters, the least precisely known is (v 2 )j/^,. Since this is a 



subtracted quantity, it can be either positive or negative. There is a usefu 



relation, first 



derived by Gremm and Kapustin using equation of motion of NRQCD 40j, relating this 
quantity with the pole mass of charm quark and the J/ip mass 8 . To our purpose this 



relation reads 



25|] 



^ m cpole 



4m cpole 



(25) 



For a reformulation of Gremm-Kapustin relation from a even lower energy effective theory of NRQCD, 
dubbed potential NRQCD, we refer the interested readers to Ref. [4l|. 
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In the analysis of the process e + e~ — > J/ipVc at B factory 25j, the charm quark pole mass 
is taken as the commonly quoted value, 1.4 GeV, consequently (v 2 )j/^ = 0.22 and (v 2 )^ = 
0.13. However, since m cpo \ e contains renormalon ambiguity, it cannot be determined better 
within an accuracy of order Aqcd- Therefore one should not be surprised that somewhat 
larger value of m cpo \ e may be occasionally reported in literature. For example, a recent QCD 
moment sum rule analysis claims m cpo ie = 1.75 ± 0.15 GeV [421 ] . For this specific value of 
charm quark pole mass, one would obtain instead (f 2 )j/^ = —0.22 ± 0.15. 

The situation seems to be rather obscure since even the sign of (v 2 )j/^ cannot be un- 
ambiguously determined. To proceed, we take a practical attitude and adopt the value 
( v2 )j/4> = 0-25 ± 0.09 from a recent Cornell potential model based analysis [43j. This pos- 
itive value seems to be most welcomed to alleviate the alarming discrepancy between the 



predicted and the measured cross section for e + e~ — > J/iJj + r\ c [25j, also seems compatible 
with a recent QCD sum rule determination of (v 2 )^ 44(. However, we leave the possi- 
bility open that this quantity may turn out to be negative after future scrutiny. Taking 
m b = M Vb /2 « 4.7 GeV, m c = 1.5 GeV, a = 1/137 and a s {m b ) = 0.22, we then find 9 

Br[^ -> Jji) J/ip] = 2.4itg x 10~ 8 . (26) 

The uncertainty is estimated by varying mb and m c in the ±100 MeV range, varying a s (fj,) 
in the —0.04 range (which corresponds to slide the scale from fi = nib to 2rrib), as well 
as taking into account the errors in measured T e + e - (of ±0.14 keV) and in (v 2 )j/^. The 
constructive interference between electromagnetic and strong amplitudes has modest effect, 
i.e., neglecting the QED contribution decreases the branching ratio by few to ten per cents. 

Our prediction for the branching ratio is at least three order-of-magnitudes smaller than 
the BFL estimate, ([T]). Despite the large uncertainties inherent in various input parameters, 
we believe our prediction captures the correct order of magnitude, 10~ 8 . It is also worth 
noting that our prediction for this exclusive decay ratio is about one thousandth of the 
inclusive 4-charm ratio, (T5]), which seems fairly reasonable. 

Experimentally J/ip can be tagged cleanly by its decay to muon pair. Multiplying (126]) 
by the branching ratios of 6% for each of the decay J/ij) — > we obtain Br[^ — > 



9 Note the actual prediction, Br ~ 10~ 8 , is much larger than the expectation based on the asymptotic 
scaling behavior Br ~ ^ »J°(m c /m(,) 8 ~ 10 -U . This can be attributed to the large prefactors in the 
right hand side of 
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J/i/jJ /if) — > 4/x] ~ (0.2 — 2.4) x 10~ 10 . The total cross section for rjb production at Tevatron 
energy is about 2.5 fib [20]. Therefore, the production cross section for this 4/i decay mode 
is about 0.05 — 0.6 fb. For the full Tevatron Run 1 data of 100 pb _1 , we then obtain between 
0.005 and 0.06 produced events. We now can safely assert that the seven 4/i events reported 
by Ref. [lj must come from sources other than rjb decay. 

Tevatron Run 2 plans to achieve an integrated luminosity of 8.5 fb -1 by 2009. Assuming 
equal a(pp — > rjb + X) at <Js = 1.96 and 1.8 TeV, we then estimate there are about 0.4-5 
produced events. Since the kinematical cuts, as well as taking into account the acceptance 
and efficiency for detecting muon, will further cut down this number, we conclude it is not 
realistic to search rjb through this decay mode in Tevatron Run 2. 

To fathom the observability of rib through this mode at LHC, we need first know the 
inclusive rjb production rate. There are rough estimates for the Xbo,2 cross sections at LHC, 
which are about 6 times larger than the corresponding cross sections at Tevatron 34]. 
Assuming the same scaling also holds for rjb, we then expect the cross section for rjb at LHC 
to be about 15 fib, subsequently the production cross section for the 4 ji events to be about 
0.3-3.6 fb. For a 300 ftT 1 data, which is expected to be accumulated in one year running 
at LHC design luminosity, the number of produced events may reach between 100 to 1000. 
The product of acceptance and efficiency for detecting J /if) decay to fi^ fT is estimated to 
be e ~ 0.1 which is perhaps a conservative estimate for LHC. Multiplying the number 
of the produced events by e 2 , we expect between 1 and 10 observed events per year. If we 
loose the constraint that J /if) must be tagged by fi + fi~ pair and also allow its reconstruction 
through e + e~ mode, we can have 4-40 observed 4-lepton events per year. 

The above analysis seems to indicate that, the chance of observing rjb at LHC through the 

4-lepton mode subsists, but critically hinges on whether the signal events can be singled out 

from the abundant background events. The most important background events may come 

from the direct double J/if> production from gg fusion 47, 4§|. From previous analysis, we 

know that all of the seven 4 ft candidate events selected by CDF based on Tevatron Run 1 
I L 

data |19[], should be regarded as this kind of background events, which seem to outnumber 
the expected signal events by several orders of magnitude. The same situation may also 
apply to LHC. It might be possible for experimentalists to judiciously choose kinematical 
cuts to significantly suppress the background events while retaining as many signal events 
as possible. 
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V. OTHER DECAY CHANNELS OF rj b 



The decay channel we have considered so far, rjb — > J/ipJ/ip, which has very clean signa- 
ture, is unfortunately very much suppressed because of its maximal violation of the helicity 
selection rule 10 . There are other two-body decay channels, e.g., rjb decays to two light mesons 
and to two charmed mesons, some of which do conserve the hadron helicity, thus may have 
much larger branching ratios. In this section, we attempt to estimate the order of magni- 
tude of the branching ratios for these processes. In addition, we also present some crude 
estimation for the rjb decays to three pseudo-scalar states. 

A. rib (j] c ) decays to two and three light mesons 

When contending with light mesons in hard exclusive processes, the most appropriate 
description of them is in terms of light-cone expansion approach. On the other hand, the 
constituent quark model, which treats the light mesons as non-relativistic bound states, is 
also frequently invoked as an alternative method for a quick order-of-magnitude estimate. 
In this sense, the preceding formulas derived for rjb — > J/ipJ/ip can be applied to describe 
the unnatural decay processes rjb — > VV, once we understand we are working with the 
constituent quark model. 

We take rjb — > (fxp as an representative. By regarding as a strangeonium, we can directly 
use ( I22p . only with some trivial changes of input parameters. We take the constituent quark 
mass m s « = 0.5 GeV. The wave function at the origin of <p, ipAO), can be extracted 

analogously from its measured dielectron width of 1.27 ± 0.04 keV [18] through (1241) . We 
take (v 2 )^ ~ 1 to reflect the fact that s is inherently not a heavy quark and <p is not truly 
a heavy quarkonium. Taking m& = 4.7 ±0.1 GeV, varying the strong coupling constant 
between a s {mb) = 0.22 and a s {2mb) = 0.18, and including the experimental uncertainty in 
r e + e -, we obtain 



The interference effect is more pronounced in this case due to larger ratio of nib to m 



This process is the only possible one for r\b decays to two ground charmonium states allowed by C and P 




(27) 



mvanance. 
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Neglecting QED contribution will decreases the branching fraction by about 20%. If this 
estimate is reliable, such a rare decay mode perhaps will never be observed experimentally. 

It is interesting also to consider the similar decay process r\ c — > 00, which has been 
measured long ago. Parallel to the preceding procedure, varying the strong coupling constant 
from a s (m c ) = 0.36 to a s (2m c ) = 0.26, taking m c = 1.5 ± 0.1 GeV, we obtain 

Bi[r] c -»• 00] = (0.3-1.5) x 10~ 5 . (28) 

The interference effect is rather modest, i.e., neglecting QED contribution only decreases 
the branching ratio by less than one per cent. This prediction is consistent with an early 



22|] 



ii 



constituent quark model based analysis using Bethe-Salpeter bound state formalism 
Note this estimate is indeed much smaller than the measured value Br exp [r] c — > 
(2.7 ± 0.9) x 10 -3 , which reflects a generic symptom in pQCD calculation, also present in 
the 7] c decays to pp and K*K* . 

As stressed several times before, some nonperturbative mechanisms must be called for to 
rescue this discrepancy. Among different proposals, a particularly attractive and predictive 



scheme has been put forward by Feldmann and Kroll (henceforth FK) [23j, by generalizing 
their influential work together with Stech on rj — rj' mixing in quark flavor basis [50j], to 
include t] c . FK models the amplitude of r/ c — > VV as the product of a factor governing 
the small admixture amplitude between r/ c and rf- 1 ' , which presumably can be inferred from 
QCD Ua(1) anomaly, times a soft vertex function parameterizing the amplitude for virtual 
rj, rj' transiting into VV. To be specific, the r\ c — > 00 in their ansatz can be described as 

where the subscript FK is used to distinguish from our earlier prediction based on the pQCD 
plus constitute quark model ansatz, Eq. (1281) . The parameter c mix depends on the specific 
flavor content of VV states but not on the initial charmonium state, whereas the charm 
mass dependence is encoded in 

9 mix (Vc) = }~ F v{l)vv (s = Ml) = ±- F v( , )vv (s = 0) ( M f_ K a ) , (30) 

JT}c J rj c \ :q c / 



in Table 1 of Ref. 



49] 



The predictions of Ref. 22[ , which employ three different forms of BS wave functions for V, are duplicated 
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where f Vc is the decay constant of r] c , the second factor indicates the on-shell coupling 
rj^VV, and the last factor parameterizes the s— dependence of this vertex function. A is 
the cutoff of typical size around 1 GeV, and the variable n = 1 corresponds to the familiar 
monopole ansatz for form factor, and n = 2 characterizes the dipole ansatz. One noteworthy 
fact is that, this simple model seems able to account for both the absolute and the relative 
strengths of partial widths for each decay channel (VV = pp, K*K* and 00), to a reasonably 
satisfactory degree. 

One curious question is to ask, whether the FK mixing mechanism can be extrapolated 
to r]b — > VV, and if so, whether the corresponding prediction to the branching ratios differs 
drastically from our pQCD-based estimate. Let us now examine this question. After some 
straightforward derivation, it is easy to find 

5£5i^4 „ U _ i^V 1/2 ( W 1 = 0.16 x 0.01", (31) 

where we have resorted to heavy quark spin symmetry for decay constants f Vc ~ jji^ and 
f Vb w / x , and consequently taken f Vc = 417 MeV and f Vb = 715 MeV through the relation 



fj/i> — y ifjj/^O). For simplicity, we have utilized the hierarchy rrib } m c ^> A in 

deriving (13TT) . If Brpx^c — > 4>4>] is identified with the measured value 2.7 x 10~ 3 , we then 
obtain BrpK^b — * 00] to be about 4 x 10~ 6 for n = 1 and 4 x 10 -8 for n = 2. It is easy to 
see that the Br FK [?7;, — > VV] scales as m^ 4( ' 1+ri ^f b " 6 , since f Vb ~ . This indicates this 

nonperturbative mechanism, also resemble the pQCD analysis in that the branching ratio is 
highly suppressed with respect to "leading twist" scaling oc 1/mf. Despite large uncertainty, 
the predicted branching ratios of rjf, — ► VV in this mixing ansatz are (much) larger than our 
earlier predictions based on "NRQCD" plus constitute quark model, (J271) . It is reasonable to 
question the validity of the above nonperturbative mechanism at rji, energy, in our opinion, 
perhaps the pQCD result in (f27j) is more trustworthy. At any rate, it will be of interest for 
future experiments to distinguish these different mechanisms, though the task of recording 
this rare decay mode in hadron machine could be extremely challenging. 

One related class of decay channels may deserve some attention. The helicity-conserving 
decay, such as r] c —>■ K»K ("||" implies longitudinally-polarized), has still not been observed 
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yet. The experimental bound for this channel is 18 
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Br exp [?? c -> K*K + c.c] < 1.28%. (32) 

Since the coupling rf'^VP is negligible, it is reasonable to assume that the aforementioned 
nonperturbative mechanism responsible for r] c — > VV, may not play a significant role for 
r] c — > VP. One may then make use of the scaling behaviors derived earlier to interconnect 
each of them. Firstly we expect that flavor SU(3) is respected in pQCD calculation of 
T] c — > VV, accordingly our prediction is Br[r/ C — > K*K ] ~ Br[r] c — > 00] ~ 10 -5 . As we 
have learned from Sec. [Til & t tree level, the maximal helicity violation in r)b — > J/ipJ/ip is 
entirely due to the charm quark transverse momentum, as is evident in ([6]). We assume the 
similar pattern occurring here, i.e., the branching ratios of r) c — > VV are also suppressed by 
(q±/m c ) 4 relative to those of the corresponding helicity-conserving channels. We then have 

4 



Br[ Vc - ITK] ~ Br[^ c - 00] ( ^p- j 6 2 su(3) ~ lO" 3 x 6 2 su(3) , (33) 



m c 

^QCD, 

where we have taken qj_ ~ Aqcd — 500 MeV. We have intentionally included a parameter 
£su(3)) to embody the extent of SU(3) flavor violation at the amplitude level. Obviously, 
this parameter should be proportional to current quark mass difference m s — md, u , divided 
by some typical hadronic scale (presumably independent of m c ). In contrast to the isospin 
violation effect, the U-spin violating process should not receive too severe suppression. If 
assuming somewhat optimistic value esu(3) ~ 0.3, we then obtain Br[?7 c — > K*K] ~ 10~ 4 , 
which is safely below the experimental bound, ( l32|) . Needless to say, it will be very beneficial 
for the ongoing high-luminosity charmonium facility such as BES III to impose more tight 
constraint on this branching ratio. 

Carrying over the same line of argument to rjb decay to VP, we obtain 

Brfr 6 - K*K] ~ Br[ Vb - 00] e 2 su(3) ~ lO" 5 x e 2 u(3) . (34) 

If we again take Cgu(3) ~ 0.1, the branching ratio is estimated to be around 10~ 6 , about one 
hundred times larger than Br[?7 b — > J/ijj J/ip}. Of course, one should not take this crude 
estimate too seriously, and a concrete pQCD calculation based on light-cone expansion 
approach, which incorporates m s and m d difference, might be illuminating. If this estimate 



12 Although this decay channel is favored by helicity selection rule, it violates the U-spin conservation [tl.!^. 
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is trustworthy, one then expects there are already about O(10 2 ) produced events in the 
Tevatron Run 1. Since K* almost exclusively decays to Ktt, one needs to select those 
resonant KKn events. However, there are practical problems about usefulness of this kind 
of decay mode in hadron collision experiments. Because of copiously produced kaons and 
pions in a typical hadron collision, huge combinatorial backgrounds might make it very 
difficult to identify the true signal. On the other hand, the prospective Super B factory, 



which runs at several T(nS) resonances with an unprecedented luminosity 51|, may produce 
an enormous amount of rjb through Ml transition from T(nS) states. With much suppressed 
backgrounds, Super B factory may offer a viable environment to detect this decay mode. 

Finally we exploit one experimental fact that some 3-body channels, i.e., rj c decays to 
three pseudo-scalar states, have exceedingly large branching ratios. To be concrete, three r\ c 
decay modes with largest branching ratios are 18] 



Bt cxp [ Vc -> KKir] = (7.0 ±1.2)%. (35) 
Br cxp [7/ C -^vrvr] = (4.9 ±1.7)%. (36) 
Br exp [ry c t/tttt] = (4.1 ±1.8)%. (37) 

Although we do not know how to calculate these processes in practice, some general pattern 
may still be identified. Subtracting off the phase space effects, one finds these three ampli- 
tudes roughly respect the SU(3) flavor symmetry. This may be taken as the sign that these 
processes could in principle be described by the pQCD scheme. These processes proceed as 
two steps. Firstly t] c annihilates to two highly virtual gluons, which then transit into four 
energetic light quarks. This is a short- distance process. Subsequently these light quarks ma- 
terialize into three pseudo-scalars via soft nonpertubative dynamics, which is a long-distance 
process. Factorization is expected to hold between these two stages due to the hard scale set 
by large c mass. Because of their large decay ratios, these processes are naturally expected 
to possess the "leading twist" scaling behaviors, Br ~ l/ m c- This scaling assumption can 
then be used to infer the ratios for rjb decays to the same pseudo-scalar states. For example, 
we may expect 

Br[r] b -> KKn] ~ Br exp [r) c KKn] (^) ~ 10" 4 . (38) 

This is by far the largest branching ratio we have found in all exclusive hadronic decay 
channels of rjt,. If this is the case, these channels will be worth pursuing. Practically, 
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the rjirir mode may not be very useful in hadron hadron collider experiments, since the 
ordinary way of reconstructing rj, which goes through the 2j decay, may suffer enormous 
contamination from combinatorial background. For similar reason, the K + K~7r° channel 
may also be difficult to detect due to presence of tt°. In contrast, the decay mode KsK^ti t 
is much more advantageous since K$ can be reconstructed cleanly via its decay to 7r + 7r~ . In 
any event, all these decay channels have promising observation potential at future Super B 
factory. 

B. 7/5 decays to two charmed mesons 

Heavy rjb mass opens the gate for it to decay into charmed particles as well as light ones. 
In fact, Maltoni and Polosa have suggested that rjb decays to two charmed mesons could be 
the most promising channels for observing rjb in Tevatron Run 2 |20j. They first perform a 
perturbative calculation for rjb decay to the inclusive ccg state, 

Br[ Vb ^ccg] = 1.5±°j%, (39) 

then interpret this value as an upper limit for the ratios of the exclusive decays to double 

u\ 

charmed mesons. MP continue to argue that, the exclusive decays to D*D (the charge 
conjugate is implicitly implied) may saturate the inclusive charmful decay ratio, and subse- 
quently estimate 1(T 3 < Bi[rj b -> D*T3 { * ] ] < 1CT 2 . 

In our opinion, this saturation assumption seems to be physically unwarranted, conse- 
quently MP's predictions for the ratios may well be an overestimate. Firstly, there seems 
no reason to expect that the parton process rjb — > ccg will be dominated by the two-body 
exclusive decays, since the g jet may readily hadronize in an independent direction, which 
will result in a multi-body decay configurations. This can be exemplified by the fact we just 
mentioned, that rj c decays to 3 pseudo-scalars have larger branching ratios than any 2-body 
decay channels. One may notice at the lowest order in a s , rjb decay to double charm mesons 
is also depicted by Fig. [U with one cc pair replaced by a qq pair. The gluon which is on 
shell in the inclusive r]b — > ccg process, now has to convert to a light quark pair with large 
invariant mass, so is highly virtual, the corresponding ratios of double charm mesons thus 
must be at least suppressed by one factor of a s and one factor of l/mf relative to (1391) . 
Taking into account the nonperturbative binding probability, which is much less than 1, will 
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further suppress the exclusive 2-body decay rates. Moreover, there seems also no strong 
reason to believe that the binding of c with q will necessarily be saturated by the ground 
state charm meson only. 

Despite the lack of an explicit pQCD calculation, one may still proceed with some physical 
consideration. The branching ratios of 775 decays to two charm mesons can depend on three 
dimensional parameters: rrib, m c and Aqcd- Since the decay rjb — > DtD conserves the 
helicity, we thus expect the corresponding branching fraction scales as 1/mf. For each 
pair of c and q to form a meson, there is a factor proportional to AqcdA^c which 
characterizes the binding probability [45!]. This is the only place where Aqcd dependence 
can enter. Therefore, we expect the following asymptotic behavior: 

Br[ Vb - D*D] ~ a] (^^) ~ 10~ 5 , (40) 

where m 2 c is inserted to make the ratio dimensionless. Notice this value is reasonably com- 
patible with the inclusive upper bound, (1391) 13 . 

We next turn to the other charmful decay channel, rjb — > D* ± Dx (rjb — * DD is forbidden 
by P invariance), which violates the helicity selection rule maximally. Assuming again the 
helicity violation is caused by the transverse momentum of light quark in D meson, the same 
as in rjb —> J/ipJ/ip, we then obtain 




where q± ~ Aqcd — 500 MeV is used. Rather small branching ratio renders this decay mode 
not so useful for detecting rjb. Nevertheless, it might be possible that a large prefactor may 
arise in the actual calculation, like what occurs in rjb — ► J/ipJ/ip (see footnote [9]) , so that 
the actual ratio may be somewhat larger than this naive estimate. In any event, a reliable 
pQCD calculation of the decay rates for rjb — > D*D will be helpful. 

Based on the previous discussions about 3-body decay of rjb, one may wish that the 
branching ratio for r] b — ► DDtx might be as large as 10~ 4 , similar to that for r]b — > KKn. 
However, this decay mode may not be as competitive as KsK ± 7r T mode, since D meson 
does not possess a clean signature comparable with K$- 

13 However, as will be reported in Ref. the decay rate of this process vanishes if D and D* are exactly de- 
generate. Therefore, the non-vanishing decay amplitude must be induced by heavy-quark-spin-symmetry 
violating effects, oc Aqcd/wc- This indicates the estimate (HPT) may well be subject to further suppression. 
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To summarize, the rjb — ► D*D channel, with an expected branching ratio ~ 10~ 5 , may 
be regarded as a valuable searching mode. The D*°D° channel (charge conjugate state 
implicitly included) may not be so useful, since D*° predominantly decays to D°7r° and D 0r y, 
where neither 7r° nor 7 can be cleanly tagged in hadron collision environment. In contrast, 
it is more advantageous to use D* + D~ mode as a trigger. D* + can be tagged through 
its decay to D°tt + , subsequently D° may be reconstructed from K~tt + , while D~~ can be 
tagged through its decay to K + tt~tt~. It is worth emphasizing that due to the proximity 
of D* + mass to the sum of masses of D° and 7r + , D* + can be cleanly identified with a 
rather narrow peak in Dtt invariant mass spectrum. Using the measured branching ratios 
Br[L>*+ -> D°n + ] « 70%, Br[D° -> R-tt + ] « 4% and Br[D" -> K + it-it-] « 10% [18], we 
estimate Br [775 — > D* + D~ — ► i^ + fr~7r + 7r~7r + 7r~] ~ 10" 8 . There are roughly 0(1) produced 
events in Tevatron Run 1, about O(10 2 ) produced events in Run 2, and about O(10 4 ) 
produced events in one year run at LHC. The statistics seems to be enough in the forthcoming 
hadron collider program for observing rjb through this decay mode, provided that the signal 
events are not swallowed by the possibly large combinatorial backgrounds. 

VI. SUMMARY AND OUTLOOK 

The major motif of this work is to clarify a controversy about whether double J /if} can be 
a useful decay mode to detect rjb in Tevatron Run 2. We have shown this process is subject to 
large kinematical suppression due to the maximal violation of the helicity selection rule. By 
an explicit pQCD calculation based on NRQCD, we predict the corresponding branching 
ratio to be only of order 10~ 8 , thus making the search for rjb through this mode rather 
unrealistic in Tevatron Run 2. Nevertheless, we anticipate that at LHC, the observation 
potential of this decay mode may not be so pessimistic, if experimentalists can find a good 
way to suppress the rather copious QCD background events. 

The large mass of rjb renders any of its exclusive decay channels in general very small. To 
provide some useful guidance for experimental search for this elusive particle, it is valuable 
to identify those decay modes with relatively large branching fractions and also with clean 
signature. To make progress along this direction, in the following we outline some issues 
which we think may deserve further studies: 

1. Stimulated by rather large QCD radiative correction to exclusive double charmonium 
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production at B factory [30l . l3lj . one may inquire how large the effect of the next- 
to-leading order QCD correction to rjb — > J/ipJ/tp is. It turns out that, at NLO in 
a s for this process, the helicity selection rule can be violated by finite charm mass, 
consequently one obtains the non-vanishing result even at the the leading order in 
velocity expansion 52] . At the amplitude level, the ratio of the radiative correction 
piece to the relativistic correction piece as considered in this work, (|20|) . is about 
— : (v 2 )j/<p ~ which implies both effects are equally important, and the more 

accurate prediction will crucially depend on their relative phase. To this end, a precise 
determination of the quantity (v 2 ) j/^, in particular its sign, would be helpful. 

2. The process rjb — > K*K+c.c, favored by the helicity selection rule but at the same time 
violating U-spin symmetry, is estimated to have a branching fraction ~ 10~ 6 -10 -7 . 
It might be worthwhile if an actual pQCD calculation which implements m s — rrid 
difference, can be carried out in the light cone expansion scheme, to compare with 
this rough estimate. Moreover, the eventual experimental sighting of r) c — > K*K in 
charmonium factory like BES III will definitely enrich our understanding toward this 
class of helicity-conserving yet flavor SU(3) violating quarkonium decay processes. 

3. The individual decay modes of r] c with largest branching ratios are i] c — > KKtt, i]7i7i 
and rfinr. Stimulated by this experimental fact, one may hope that the 3-body decay 
channels of rjb, such as KsK ± ir T , with an expected branching ratio of order 10~ 4 , 
might be a potentially useful searching mode for r)b in current and forthcoming hadron 
collider programs, if it can survive from the copious combinatorial background events. 
This mode will definitely have promising potential to be observed in the prospective 
Super B factory. 

4. Another helicity-conserving decay process, r\ h — > D*D + c.c, with an expected branch- 
ing ratio of order 10~ 5 , may also be well worth searching for experimentally. It is 
also profitable to carry out a concrete calculation of this exclusive double charm decay 
process from pQCD scheme, but one may be obliged to incorporate the heavy-quark- 



spin-symmetry violating effect 46]. 



5. In this work we haven't discussed the possibility of observing rjb through its decay 
to baryon pair, such as r\b — > pp. If the corresponding ratio is not too small, this 
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is potentially a good searching mode thanks to relatively fewer baryonic background 
events in hadron collision experiments 14 . Experimentally rj c — > pp is observed to 



have a branching ratio of order 10 3 18]. However, one should be aware that this 



process also violates the helicity selection rule, and the available pQCD prediction, 
when taking into account the nonzero light quark mass but still in the collinear quark 



configuration, is still far smaller than the measured value 53J, therefore some nonper 



turbative mechanism needs to be called for to explain this discrepancy [23]. Because 
of rather heavy mass of rjf,, one may hope pQCD framework should provide a reliable 
prediction for rjb — > pp. It will be valuable if a more systematic calculation can be 
carried out. 

Note added. While the revised version of the manuscript is to be submitted, there 



appears an eprint in arXiv by Babar collaboration [54|, which claims the first unambiguous 
discovery of % through the hindered Ml transition process T(3S) — > 77^7. It is interesting 
to note the rather precisely measured 775 mass, 9388.9±H(stat) ± 2.7(syst) MeV, seems not 
compatible with the predictions from most potential models as well as the weakly-coupled 



potential NRQCD, instead consistent with the lattice QCD prediction within error 
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